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Cosmic rays are hampered by the Moon and a deficit in its direction is expected (the so-called 
Moon shadow). The Moon shadow is an important tool to determine the performance of an air 
shower array. Indeed, the westward displacement of the shadow center, due to the bending effect of 
the geomagnetic field on the propagation of cosmic rays, allows the setting of the absolute rigidity 
scale of the primary particles inducing the showers recorded by the detector. In addition, the shape 
of the shadow permits to determine the detector point spread function, while the position of the 
deficit at high energies allows the evaluation of its absolute pointing accuracy. 

In this paper we present the observation of the cosmic ray Moon shadowing effect carried out 
by the ARGO-YBJ experiment in the multi-TeV energy region with high statistical significance 
(55 standard deviations). By means of an accurate Monte Carlo simulation of the cosmic rays 
propagation in the Earth-Moon system, we have studied separately the effect of the geomagnetic 



field and of the detector point spread function on tiie observed sliadow. Tlie angular resolution as 
a function of the particle multiplicity and the pointing accuracy have been obtained. The primary 
energy of detected showers has been estimated by measuring the westward displacement as a function 
of the particle multiplicity, thus calibrating the relation between shower size and cosmic ray energy. 
The stability of the detector on a monthly basis has been checked by monitoring the position and the 
deficit of the Moon shadow. Finally, we have studied with high statistical accuracy the shadowing 
effect in the "day" /"night" time looking for possible effect induced by the solar wind. 

PACS numbers: 96.50.S-;96.50.sd;96.50.Bh 



I. INTRODUCTION 

Cosmic rays (CRs) blocked in their way to the Earth 
by the Moon generate a deficit in its direction usually 
mentioned as 'Moon shadow '. The analysis of the Moon 
shadow observed by an air shower array may provide 
unique information on its performance. At high ener- 
gies, the Moon shadow would be observed by an ideal 
detector as a 0.52° wide circular deficit of events, cen- 
tered on the Moon position-'^. The actual shape of the 
deficit as reconstructed by the detector allows the deter- 
mination of the angular resolution while the position of 
the deficit allows the evaluation of the absolute pointing 
accuracy. In addition, charged particles are deflected by 
the geo magnetic field (GMF) by an angle depending on 
the energy. As a consequence, the observation of the dis- 
placement of the Moon shadow at low rigidities can be 
used to determine the relation between the shower size 
and the primary energy. 

The same shadowing effect can be observed in the di- 
rection of the Sun but the interpretation of this phe- 
nomenology is less straightforward. In fact, the displace- 
ment of the shadow from the apparent position of the 
Sun could be explained by the joint effects of the GMF 
and of the solar and interplanetary magnetic fields (SMF 
and IMF, respectively), whose configuration considerably 
changes with the phases of the solar activity cycle [l|. 
In this regard, understanding the Moon shadow phe- 
nomenology is a useful tool to investigate the GMF fea- 
tures needed to disentangle the effects of different mag- 
netic fields on the Sun shadow and to perform a measure- 
ment of the IMF during a minimum of the solar activity 

i- 

Finally, the Moon shadow can be exploited to measure 
the antiproton content in the primary CRs. In fact, act- 
ing the Earth-Moon system as a magnetic spectrometer, 
paths of primary antiprotons are deflected in the oppo- 
site direction with respect to those of the protons in their 
way to the Earth. This effect has been used to set limits 
on the antiproton flux at TeV energies not yet accessible 
to balloon or satellite experiments [3|-|6[. 

In this paper we present the observation of the cosmic 
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^ Actually, the width of the Moon disc ranges from 0.50° to 0.58° 
depending on its distance from the Earth. 



ray Moon shadowing effect carried out by the ARGO- 
YBJ experiment during the period from July 2006 to 
November 2010. We report on the angular resolution, the 
pointing accuracy and the rigidity scale calibration of the 
detector in the multi-TeV energy region. The results are 
compared with the predictions of a detailed simulation 
of cosmic ray propagation in the Earth- Moon system. 

The paper is organized as follows. In Sec. II the 
ARGO-YBJ detector is described and the event recon- 
struction sketched out. In Sec. Ill the data analysis per- 
formed with two different background estimation tech- 
niques is outlined. The results of a Monte Carlo simu- 
lation of the cosmic ray propagation in the Earth-Moon 
system are presented in Sec. IV. The measurement of the 
pointing accuracy and of the angular resolution as well as 
the evaluation of the absolute rigidity scale are discussed 
in Sec. V. A high statistics study of the day-night effect 
is also reported in Sec. V. A summary of the obtained 
results is given in Sec. VI. 



II. THE ARGO-YBJ EXPERIMENT 
A. The detector 

The ARGO-YBJ experiment, located at the YangBa- 
Jing Cosmic Ray Laboratory (Tibet, P.R. China, 4300 m 
a.s.l., 606 g/cm^), is currently the only air shower array 
exploiting the full coverage approach at very high alti- 
tude, with the aim of studying the cosmic radiation at 
an energy threshold of a few hundred GeV. 

The detector is constituted by a central carpet ~74x 
78 m^ , made of a single layer of resistive plate chambers 
(RPCs) with ^93% of active area, enclosed by a guard 
ring partially instrumented (~20%) up to ^100x110 m^. 
The apparatus has a modular structure, the basic data 
acquisition element being a cluster (5.7x7.6 m^), made of 
12 RPCs (2.85x1.23 m^ each). Each chamber is read by 
80 external strips of 6.75x61.80 cm^ (the spatial pixels), 
logically organized in 10 independent pads of 55.6x61.8 
cm^ which represent the time pixels of the detector Q . 
The readout of 18360 pads and 146880 strips is the ex- 
perimental output of the detector. The relation between 
strip and pad multiplicity has been measured and found 
in fine agreement with the Monte Carlo prediction 7] . In 
addition, in order to extend the dynamical range up to 
PeV energies, each chamber is equipped with two large 
size pads (139x 123 cm^) to collect the total charge devel- 
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final configuration of 153 clusters since November 2007 
with the trigger condition Ntrig = 20 and a duty cycle 
>85%. The trigger rate is ~3.5 kHz with a dead time of 
4%. 

In the present study the data recorded by the digital 
readout have been analyzed to measure the Moon shadow 
effect induced by low-energy primaries. 



B. Event reconstruction and data selection 



FIG. 1: Layout of the ARGO-YBJ experiment (see text for a 
detailed description of the detector). 



oped by the particles hitting the detector !3|. The RPCs 
are operated in streamer mode by using a gas mixture (Ar 
15%, Isobutane 10%, TetraFluoroEthane 75%) for high 
altitude operation Q- The high voltage settled at 7.2 kV 
ensures an overall efficiency of about 96% [l(| . The cen- 
tral carpet contains 130 clusters (hereafter ARGO-130) 
and the full detector is composed of 153 clusters for a 
total active surface of ^6700 m^ (Fig. [1]). The total 
instrumented area is ^11000 m". The information on 
strip multiplicity and the arrival times recorded by each 
pad are received by a local station devoted to manage 
the data of each cluster. A central station collects the 
information of all the local stations. The time of each 
fired pad in a window of 2 /xs around the trigger time 
and its location are used to reconstruct the position of 
the shower core and the arrival direction of the primary 
particle. In order to perform the time calibration of the 
18360 pads, a software method has been developed pl| . 
To check the stability of the apparatus a control system 
(DCS) monitors continuously the current of each RPC, 
the gas mixture composition, the high voltage distribu- 
tion as well as the environment conditions (temperature, 
atmospheric pressure, humidity). A simple, yet powerful, 
electronic logic has been implemented to build an inclu- 
sive trigger. This logic is based on a time correlation 
between the pad signals depending on their relative dis- 
tance. In this way, all the shower events giving a number 
of fired pads Npad > ^trig in the central carpet in a time 
window of 420 ns generate the trigger. This trigger can 
work with high efficiency down to Ntrig = 20, keeping 
negligible the rate of random coincidences. 

Because of the small pixel size, the detector is able 
to record events with a particle density exceeding 0.003 
particles m~^, keeping good linearity up to a core density 
of about 15 particles m~^. This high granularity allows 
a complete and detailed three-dimensional reconstruction 
of the front of air showers at an energy threshold of a few 
hundred GeV. Showers induced by high energy primaries 
(> 100 TeV) are also imaged by the analog readout of 
the large size pads ^8j . 

The whole system, in smooth data taking since July 
2006 with ARGO-130, is in stable data taking with the 



The reconstruction of the shower parameters is carried 
out through the following steps. 

At first, a plane surface is analytically fitted (with 
weights equal to 1) to the shower front. This proce- 
dure is repeated up to 5 times, each iteration rejecting 
hits whose arrival time is farther than 2 standard de- 
viations from the mean of the distribution of the time 
residuals from the fitted plane surface. This iterative 
procedure is able to reject definitively from the recon- 
struction the time values belonging to the non-Gaussian 
tails of the arrival time distributions ^2]] . After this first 
step the problem is reduced to the nearly- vertical case by 
means of a projection which makes the fit plane overlap- 
ping the detector plane. Thereafter, the core position, 
i.e. the point where the shower axis intersects the detec- 
tion plane, is obtained fitting the lateral density distribu- 
tion of the secondary particles to a modified Nishimura- 
Kamata-Greisen (NKG) function. The fit procedure is 
carried out via the maximum likelihood method [l3| . Fi- 
nally, the core position is assumed to be the apex of a 
conical surface to be fitted to the shower front. The slope 
of such a conical correction is fixed to a = 0.03 ns/m [IJ] . 

The capability of reconstructing the primary arrival 
direction can be further enhanced by applying robust 
statistical methods in the analysis of the shower front, 
conveniently weighting the contribution of the most de- 
layed particles. In detail, we first fit a conical surface to 
the shower image, by minimizing the sum of the squares 
of the time residuals. At this stage, all the particles 
hitting the detector have the same weight Wi=l. Af- 
ter computing the RMS of the time residual distribu- 
tion with respect to such a conical surface, we set K = 
2.5-RMS as a 'scale parameter' and perform the min- 
imization of the square of the time residuals weighted 
sum, where Wi—l if the particle is onward the shower 
front, 
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tj )/K) otherwise. The function 



f{x) is a common Tukey biweight function [14]. The 
fit procedure is iterated, every time refreshing the scale 
parameter, until the last reconstructed direction differs 
from the previous one for less than 0.1°. 

The analysis reported in this paper refers to events se- 
lected according to the following criteria: (1) more than 
20 strips Nstrip should be fired on the ARGO-130 car- 
pet; (2) the zenith angle of the shower arrival direction 
should be less than 50°; (3) the reconstructed core posi- 
tion should be inside an area 150 x 150 m^ centered on the 
detector. After these selections the number of events an- 



alyzed is about 2.5 x 10^^ (about 10^ inside a 10°x 10° an- 
gular region centered on the Moon position). According 
to simulations, the median energy of the selected protons 
is E50 «1.8 TeV (mode energy «0.7 TeV). 



III. MOON SHADOW ANALYSIS 

For the analysis of the shadowing effect three different 
sky maps in celestial coordinates (right ascension R.A. 
and declination DEC.) are built: the event and back- 
ground maps with 0.1° x 0.1° bin size, necessary to deter- 
mine the deficit shape, and the significance map used to 
estimate the statistical significance of the observation. 

The event map, of size 10° x 10°, centered on the Moon 
location, is filled with the detected events. Cosmic rays 
blocked by the Moon have to be as many as the back- 
ground events lying within a region as large as the Moon 
disc. A suitable background estimation is therefore a 
crucial point of the analysis. The background has been 
evaluated with both the time-swapping [15| and the equi- 
zenith angle [16j methods in order to investigate possible 
systematic uncertainties in the background calculation. 

In the time-swapping method, N "fake " events are gen- 
erated for each detected one, every time replacing the 
measured arrival time with a new one. Such a random 
time spans over a 3 hours wide buffer of recorded data, to 
minimize the systematic effects induced by environment 
variations (i.e. temperature and atmospheric pressure). 
Changing the time, each fake event maintains the same 
declination of the original one, but has a different right 
ascension. In this way, a new sky map (the background 
map) is built. If the number of fake events generated 
for each event is N, the fluctuations of the background 
estimation are reduced of a factor ^ \N with respect 
to those of the signal. In this analysis we set N = 10. 
The strong point of the time-swapping technique is that 
it takes into account only the sky region where the Moon 
actually passes through, though a few tens of minutes 
before or later. On the other hand, also the time of all 
the events is swapped, then the obtained background at 
the Moon position is slightly underestimated and thus 
the signal is underestimated. This underestimate ranges 
from ~4% to 10%, increasing with the angular resolution, 
hence depending on the event multiplicity. The observed 
event rate is then corrected using the appropriate factor 

With the equi-zenith angle method the number of cos- 
mic rays recorded in the off-source cells with the same 
size, at the same zenith angle and in the same time in- 
tervals as the on-source cell is averaged. The method 
is able to eliminate various spurious effects caused by in- 
strumental and environmental variations, such as changes 
in pressure and temperature that are hard to control and 
tend to introduce systematic errors in the measurement. 
The equi-zenith background estimation is achieved in the 
reference frame of the experiment, i.e. using the local 
coordinates zenith and azimuth. The Moon position is 



computed every minute and 6 off-source bins are sym- 
metrically aligned on both sides of the on-source field, at 
the same zenith angle. The nearest off-source bins are set 
at an azimuth distance 5° from the on-source bin. The 
other off-source bins are located every 5° from them. The 
average of the event densities inside these bins is taken 
to be the background. The equi-zenith technique uses 
only showers detected at the same time of the on-source 
events and permits to take into account every minimal 
sudden environment change. Nonetheless, its efficiency 
relies on the assumption that the events triggering the 
detector are uniformly distributed in azimuth, which is 
true only at the first order. As a matter of fact, different 
factors can induce a modulation in the event distribution. 
The GMF, for example, induces a modulation as large as 
^1% for low energy showers [i3;l3j making necessary a 
proper correction to the background estimation. 

The significance map is obtained from the event and 
background maps after applying the following smoothing 
procedure to take into account the angular resolution of 
the detector. The bins of the maps are 'integrated' over 
a circular area of radius tp, i.e. every bin is filled with the 
content of all the surrounding bins whose center is closer 
than tp from its center. The value of V' is related to the 
angular resolution of the detector, and corresponds to the 
radius of the observational window that maximizes the 
signal to background ratio, which in turn depends on the 
event multiplicity: when the point spread function (PSF) 
is a Gaussian with RMS cr, then ij: = 1.58-(T and contains 
^72% of the events. The optimal size of the observational 
window as a function of the event multiplicity is obtained 
from the analysis of the event map and compared with 
the results of a Monte Carlo simulation (Sec. V, B). 

After such a smoothing procedure, an integrated 
"source map" is obtained by subtracting the integrated 
background map content from that of the integrated 
event map. The deficit significance of each bin of the 
source map with respect to the content of the correspond- 
ing background map bin is computed according to Li and 
Ma [19], providing the "significance map". This map is 
used to estimate the statistical significance of the obser- 
vation. 

A detailed study of the two background calculation 
methods in the same sky region has shown that on aver- 
age they give significances of the deficit consistent within 
about 1 standard deviation, corresponding to a few per 
cent of uncertainty on the number of events in the ob- 
served Moon shadow signal. 

In the following the results obtained with the equi- 
zenith method are shown. 

In Fig. [2] the significance map of the Moon region 
observed with data recorded in the period July 2006 - 
November 2009 (about 3200 hours on-source) is shown for 
events with fired strips ^strip > 100. The opening angle 
ip used in the smoothing procedure is 1°. The statistical 
significance of the maximum deficit is about 55 standard 
deviations. The ARGO-YBJ experiment is observing the 
Moon shadow with a significance of about 9 standard 




FIG. 2: Significance map of the Moon region for events with 
^ strip >100, observed by the ARGO-YBJ experiment in the 
period July 2006 - November 2009 (about 3200 hours on- 
source in total). The coordinates are R.A. a and DEC. 5 
centered on the Moon position (om, 5m)- The color scale 
gives the statistical significance in terms of standard devia- 
tions. 



IV. MONTE CARLO SIMULATION 

A detailed Monte Carlo simulation has been performed 
in order to propagate the CRs in the Earth-Moon sys- 
tem [201 • The air shower development in the atmosphere 
has been generated with the CORSIKA v. 6.500 code 
|2l| . The electromagnetic interactions are described by 
the EGS4 package while the hadronic interactions above 
80 GeV are reproduced by the QGSJET-II.03 and the 
SYBILL models. The low-energy hadronic interactions 
are described by the FLUKA package. Cosmic ray spec- 
tra have been simulated in the energy range from 10 
GeV to 1 PeV following the relative normalization given 
in 22]. About 10^ showers have been sampled in the 
zenith angle interval 0° - 60°. The secondary particles 
have been propagated down to cut-off energies of 1 MeV 
(electromagnetic component) and 100 MeV (muons and 
hadrons). The experimental conditions (trigger logic, 
time resolution, electronic noises, relation between strip 
and pad multiplicity, etc.) have been taken into account 
via a GEANT4-based code [2j]. The core positions have 
been randomly sampled in an energy-dependent area 
large up to 2-10^ x 2- 10'^ m^, centered on the detec- 
tor. Simulated events have been generated in the same 
format used for the experimental data and analyzed with 
the same reconstruction code. 



A. The geomagnetic model 
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FIG. 3: Deficit of CRs around the Moon position projected 
along the R.A. direction. Showers with N strip >100 recorded 
in the period July 2006 - November 2009 are shown. 



To properly describe the Moon shadowing effect, the 
magnetic field from the Moon to the Earth must be taken 
into account as much accurately as possible. Since the 
contribution due to the Moon itself is negligible, the to- 
tal field turns out to be the superposition of the IMF, 
due to the solar wind, and the GMF. The latter is by far 
the most intense acting upon the particles propagating 
in the relatively narrow region between the Moon and 
the Earth. Therefore, the observed deviation of the CR 
trajectories depends mainly on the experimental site po- 
sition relative to the GMF. 

It has been already noticed that if a primary cosmic ray 
(energy E, charge Z) traversing the GMF is observed by 
a detector placed at YangBaJing, its trajectory is bent 
along the East- West direction, whereas no deviation is 
expected along the North- South one [l|. To a first ap- 
proximation, the amount of the East- West shift can be 
written as: 



'1.58° 



deviations per month. 

As can be noticed from Fig. [3l the Moon shadow turns 
out to be a deep in the smooth CR signal observed by 
ARGO-YBJ, even without subtracting the background 
contribution. The background events are not uniformly 
distributed around the Moon, because of the non uniform 
exposure of the map bins to CR radiation. 



£;[TeV] 



(1) 



The sign is set according to the usual way to represent 
the East- West projection of the Moon maps (see Figs. [H 
and[n]). Equation ^ can be easily derived by assuming 
that the GMF is due to a pure static dipole lying in the 
center of the Earth (see Appendix). As shown below, 
Eq. ([T]) is valid for nearly vertical primaries with energy 
greater than a few TeV. To perform an evaluation of the 
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FIG. 4: Deviation induced by the GMF on protons (blue 
points) and Helium nuclei (red points) . Each point refers to a 
simulated primary. The analytical trends obtained from Eq. 
((T]) are also shown as dashed (protons) and dot-dashed (He) 
lines. 



bending effect, it is necessary to adopt a model of the 
magnetic field in the Earth-Moon system. Such a model 
provides an estimation of the coefficients of the magnetic 
field expansion in spherical harmonics. The s imp lest one 
is the so-called virtual dipole model (VDM) [23 ■ A bet- 
ter choice is the Tsyganenko-IGRF model (hereafter T- 
fGRF) [2^1, which takes into account both internal and 
external magnetospheric sources by using data available 
from spacecraft missions. We compared the effect on 
the particle trajectories of VDM and T-IGRF, finding in 
both cases non negligible differences with respect to Eq. 
([1]), which underestimates the deviation up to 10 - 15%, 
mostly for low energy primaries. Among the two models 
themselves, we observed discrepancies up to ^ 10%, cor- 
responding to 0.4° - 0.7° for sub-TeV primary energies, 
mainly due to the description of the field intensity nearby 
the Earth surface. Since the T-IGRF model takes into 
account more factors, we will refer to it hereafter. 

In Fig. |4]the actual East- West displacement obtained 
applying the T-IGRF model to the propagation of both 
protons (Z=l) and Helium nuclei (Z=2) can be appreci- 
ated. The points represent the deflection undergone by a 
nucleus propagating from the Moon to the YangBaJing 
geographical site according to the following simulation 
procedure: (1) the primary energy is sampled according 
to the spectra quoted in j22]; (2) the arrival direction 
is sampled from an isotropic distribution; (3) the events 
are spread uniformly during year 2008. In this figure 
the lines reproduce the deviation expected from Eq. ([1]) 
for both protons and Helium nuclei. The analytical ap- 
proach clearly underestimates the East- West deviation, 
in particular for sub-TeV events. Figure [5] shows the dif- 
ference of the T-IGRF-induced deviation with respect to 
the leading term 1.58°Z/_E[TeV]. From the plot (a) is ev- 
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FIG. 5: Panel (a) shows the residual displacement with re- 
spect to the analytical expectation (Eq. ([l])) along the East- 
West direction as a function of the primary energy. Panel (b) 
shows the residual displacement along the North-South direc- 
tion. The deviation is calculated by applying the T-IGRF 
model (see text). The color scale represents the number of 
showers lying on the single pixel of the figure. 



ident that the T-IGRF model predicts a deviation along 
the East- West direction greater than the one expected 
from Eq. ([T]). Altough this effect is negligible for ener- 
gies E>10 TeV, at lower energies E<1 TeV the difference 
can reach 1° or more. The plot (b) shows the difference 
along the North-South direction. Notice that unlike what 
the analytical approach would suggest, the North-South 
deviation of a primary can be non-null, being zero only 
on average. 



B. Moon shadow simulation 

By following the procedure described above, we can ob- 
tain the Moon shadow maps represented in Fig. [51 where 
the effect of folding the detector PSF and the GMF is in- 
vestigated. After the simulation, only events satisfying 
the selection criteria discussed in Sec. II, B have been 
taken into account. 

In the upper left plot the Moon disc as it would be 
observed by an ideal detector without any effect induced 
by the GMF is shown. In the upper right plot the ef- 
fect of the GMF on the ideal detector is displayed. The 
showers do not gather any more in the Moon disc. Along 
the R.A. direction (also East- West hereafter), they all 
suffer a 'negative' deviation (what we call 'westward'), 
inversely proportional to the energy. The long tail of the 
left part of the map is due to the lowest energy CRs (sub- 
TeV showers) which are more deviated. Along the DEC. 
direction (also North-South hereafter) , a significant devi- 
ation is suffered only by the least energetic primaries, all 
others propagating imperturbed. In the lower plots the 
effect of the detector PSF is taken into account, without 
and with the GMF. As it can be seen from the bottom left 
plot, the detector PSF only smears out the signal, leaving 
intact the circular symmetry, as expected. The combined 
effect of the GMF and the detector PSF is shown in the 
bottom right plot. The contribution of different cosmic 
ray primaries (protons, Helium and CNO group) to the 
Moon shadow deficit is shown in Fig. [T] Events con- 
tained in an angular band parallel to the East- West axis 
and centered on the observed Moon position, compatible 
with the multiplicity-dependent angular resolution, are 
used. According to Fig. HI the displacement of Helium- 
induced showers is expected to be greater than that of 
showers generated by proton primaries. This result is 
not evident in Fig. [71 Indeed, the analysis criteria based 
on the event multiplicity (the experimental observable) 
select the rigidity distributions shown in Fig. [51 The 
Helium rigidity spectrum exhibits a mode higher than 
that of the proton rigidity spectrum, resulting in a lower 
displacement. 



C. Role of the detector Point Spread Function 

The effects of the detector PSF and of the GMF can 
be studied separately in the East- West and North-South 




FIG. 6: The effect of folding different contributions to the 
Moon signal. Upper part of the figure: Moon as it would 
be observed by an ideal detector without GMF (left plot). 
Effect of the GMF on an ideal detector (right plot). Lower 
part: effect of the detector PSF without and with the GMF 
(left and right plot, respectively). Only the showers satisfying 
the selection criteria in Sec. II, B are shown. The color scale 
represents the number of showers lying on the single pixel of 
the figure. 
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FIG. 7: Simulated deficit counts around the Moon projected 
to the East- West axis for N strip >100. The contribution of 
different primaries to the Moon shadow deficit can be appre- 
ciated. The CNO component has been multiplied by a factor 



projections. As already noticed, if we consider the mag- 
netic deviation but not the smearing due to the angular 
resolution of the detector, the symmetry of the signal is 
broken only along the East- West direction (see Fig. [51 
upper right map). Furthermore, the North-South devia- 
tion is less than Z- 0.l7E[TeV] for 95% of CRs, making 
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FIG. 8: Rigidity distributions for events induced by differ- 
ent nuclei satisfying tlie selection criteria in Sec. II, B. The 
calculations refer to the showers of Fig. [T] 



us confident that along this direction the signal is mostly 
afl^ected by the angular resolution, which can be then de- 
termined. 

The angular width of the Moon (about half a degree) 
contributes to the spread of the signal, therefore we must 
disentangle this effect in measuring the detector angular 
resolution. Assuming a Gaussian PSF with variance a^, 
the width of the observed signal results: 



RMS = a0\ 1 + 



2^ 



(2) 



where r^ is the Moon radius. The contribution of the 
Moon size to the RMS is dominant when ag is low, i.e. 
at high particle multiplicities. For instance, the difference 
between RMS and ae is 20% if ag = 0.2°, less than 5% if 
ag >0.4°, and only 1.7% Hag = 0.7°. In Fig. [9]the effect 
of the detector angular resolution along the East- West 
projection of the Moon shadow deficit is shown. Such 
an effect determines not only the smearing, but also a 
further displacement of the signal peak due to the folding 
with the asymmetrical deflection induced by the GMF. 
The West tail of the shifted signal, indeed, has a larger 
weight compared with the sharp East edge and tends to 
pull the signal in its direction. 



V. RESULTS AND DISCUSSIONS 

A. The shape of the Moon shadow 

To get information on the detector performance the 
experimental shape of the Moon shadow for different 
shower multiplicities has been compared with the results 
of the Monte Carlo simulation of CR propagation in the 
Earth-Moon system. The analysis is carried out by us- 
ing the "source" sky maps built subtracting the back- 



2000 



1500 



1000 



500 




FIG. 9: Effect of the detector PSF along the East- West di- 
rection for proton-induced showers. The continuous black 
line represents the Moon shadow deformed by the GMF as 
it would appear to an ideal detector. The segmented red line 
takes into account also the effect of the detector PSF: the 
diplacement of the signal peak results enhanced. 



ground maps to the event ones. The deficit counts ob- 
served around the Moon projected on the East- West axis 
are shown in Fig. [TU] for 4 multiplicity bands compared 
to Monte Carlo expectations. We used the events con- 
tained in an angular band parallel to the East- West axis 
and centered on the observed Moon position. The widths 
of these bands are chosen on the basis of the Monte Carlo 
simulation so that the shadow deficit is maximized. They 
turn out to be proportional to the Nstrip-dependent an- 
gular resolution: ±2.9° in 40< Nstnp < 60, ±2.6° in 
60< Nstrip < 100, ±2.1° in 100< Nstrip < 200, ±1.6° in 
200< Nstrzp < 500. 

As an expected effect of the GMF, the profile of the 
shadow is broadened and the peak positions shifted west- 
ward as the multiplicity (i.e., the cosmic ray primary en- 
ergy) decreases. The data are in good agreement with the 
expectations from the Monte Carlo simulation discussed 
in the previous section. 

The deficit counts measured around the Moon pro- 
jected along the North-South axis for the same multi- 
plicity bins are compared to Monte Carlo expectations 
in Fig. [TTJ As expected the shape of the deficit, not 
affected by the GMF along this direction, is symmetric 
and provides information about the detector PSF and 
the angular resolution. 



B. Pointing accuracy and angular resolution 

The pointing accuracy of the detector can be deter- 
mined by observing the position of the Moon shadow 
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FIG. 10; Deficit counts measured around the Moon projected 
along the East- West axis for different multiplicity bins (black 
circles) compared to Monte Carlo expectations (red squares). 
Events contained in an angular band parallel to the East- West 
axis and centered on the observed Moon position, propor- 
tional to the multiplicity-dependent angular resolution, are 
used (see text). 
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FIG. 11: Deficit counts measured around the Moon pro- 
jected along the North-South axis for different multiplicity 
bins (black circles) compared to Monte Carlo expectations 
(red squares). Events contained in an angular band parallel 
to the North-South axis and centered on the observed Moon 
position, proportional to the multiplicity-dependent angular 
resolution, are used (see text). 
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FIG. 12: Significance map of the Moon shadow region ob- 
served by the ARGO-YBJ detector in 3200 hours on-source 
for events with Nst^jp >1000. The color scale gives the sta- 
tistical significance in terms of standard deviations. 
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FIG. 13: Measured angular resolution of the ARGO-YBJ 
detector (upward black triangles) compared to expectations 
from Monte Carlo simulation (downward red triangles) as a 
function of the particle multiplicity. The multiplicity bins are 
shown by the horizontal bars. 



cast by high energy CRs whose trajectory is negligibly 
aflFected by the GMF. The deflection of 30 TeV pro- 
tons is about 0.05° (see Eq. ([T])). For heavier nuclei 
this deflection increases but as the composition of cos- 
mic rays in this energy range is dominated by the light 
component (nuclei heavier than CNO contribute to the 
observed rate for less than 3% in the whole strip mul- 
tiplicity range [23|), we expect only a negligible effect 
on the blurring of the Moon shadow from heavy ions. 
The observed position of the Moon shadow cast by high 



energy (Nstr^p >1000, Ef^ - 30 TeV) primary CRs is 
centered in the East- West direction, as shown in Fig. 
[T2l On the other hand, according to the simulation of 
CR propagation in the Earth-Moon system, no displace- 
ment along the North-South direction at any energy is 
expected at the YangBaJing latitude. This analysis sug- 
gests that there is a residual systematic shift towards 
North independent of the multiplicity. 

The PSF of the detector, studied in the North-South 
projection not affected by the GMF (see Sec. IV, C), 
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is Gaussian for Nstrip >200, while for lower multiplici- 
ties is better described for both Monte Carlo and data 
with a linear combination of two Gaussian functions. 
The second Gaussian contributes for about 20%. For 
these events the angular resolution is calculated as the 
weighted sum of the (Tg of each Gaussian. In Fig. [13] 
the angular resolution measured along the North-South 
direction is compared to Monte Carlo predictions as a 
function of the particle multiplicity, i.e. the number of 
fired strips ISigtrip on ARGO-130. The values are in fair 
agreement showing that the ARGO-YBJ experiment is 
able to reconstruct events starting from only 20 particles 
spread on an area ^6000 m^ large with an angular res- 
olution better than 1.6°. The effect of the finite angular 
width of the Moon on the angular resolution, ruled by 
Eq. ([2]), has been taken into account. 

This measured angular resolution refers to cosmic ray- 
induced air showers. The same Monte Carlo simula- 
tion predicts an angular resolution for 7-induced showers 
smaller by ~30-40%, depending on multiplicity, due to 
the better defined time profile of the showers. 



C. Absolute rigidity scale calibration 

In order to calibrate the absolute rigidity scale of CRs 
observed by the ARGO-YBJ detector we can use the 
GMF as a magnetic spectrometer. In fact, the westward 
displacement of CRs by an angle inversely proportional 
to their energy (Eq. ([T])) provides a direct check of the 
relation between the shower size and the primary energy. 
In Fig. [14] the displacements of the Moon shadow in both 
North-South (upper plot) and East- West (lower plot) di- 
rections as a function of the particle multiplicity, i.e. the 
number of fired strips N^trip on ARGO-130, are shown. 
The rigidity scale refers to the rigidity (TeV/Z) associ- 
ated to the median energy in each multiplicity bin. 

The same Monte Carlo simulation predicts that at 
fixed multiplicity the median energy for 7-induced show- 
ers is smaller by sa30% on average. 

The observed shift is compared to the results of the 
Monte Carlo simulation of CR propagation in the Earth- 
Moon system. A shift of (0.19±0.02)° towards North 
can be observed. This displacement is independent of 
the multiplicity. Many tests on the absolute position of 
the detector, on the geometry of the experimental setup, 
on the time calibration and on the software for recon- 
struction have been carried out. The most important 
contribution to the systematics is likely due to a residual 
effect not completely corrected by the time calibration 
procedure. Further studies are under way. 

Concerning the East- West direction, the good agree- 
ment between data and simulation allows the attribution 
of this displacement to the combined effect of the detec- 
tor PSF and the GMF. Therefore, the rigidity scale can 
be fixed in the multiplicity range 20-2000 particles, where 
the Moon shadow is moving under the bending effect of 
the GMF. The Monte Carlo results are fitted by the func- 
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FIG. 14: Measured displacements of the Moon shadow (upper 
plot: North- South, lower plot: East- West) as a function of 
multiplicity (black squares). The data are compared to Monte 
Carlo expectations (red circles). In the upper plot the solid 
line is fitted to the data. In the lower plot the solid curve 
is fitted to the Monte Carlo events and the dashed curves 
show the ±10% deviation from the solid one. The rigidity 
scale refers to the rigidity (TeV/Z) associated to the median 
energy in each multiplicity bin (shown by the horizontal bars). 



tion Aa = k.(N. 



strip ) 



with K 



-10.17 and A 



-0.63, 



shown by the solid curve in Fig. [T3] To estimate the 
possible shift in particle multiplicity between data and 
simulation, as shown by the dotted curves in Fig. [T31 the 
experimental data are fitted by the same function but 
with a multiplicity shift term: 



10.17[(l-Ai?„)iV, 



strip] 



-0.63 



(3) 



as described in 28] . The parameter Ai?„ is the multiplic- 
ity shift ratio, resulting in Ai?„ = (-1-4 ± 7)%. Finally, 
the conversion from A_R„ to the energy shift ratio /S.Re is 
performed. To determine the relationship between Ai?„ 
and ARe, and to check that this method is sensitive 
to energy, six Monte Carlo event samples in which the 
energy of the primary particles is systematically shifted 
event by event in the Moon shadow simulation are calcu- 
lated [2^. These six ARe samples correspond to ±20%, 
±15% and ±8%. Finally, by assuming a linear depen- 
dence, the relation Ai?„= (-0.91 ± 0.16)-ARe is ob- 
tained. Hence, the systematic uncertainty in the abso- 
lute rigidity scale ARe is estimated to be (±5 ± 8)%, 
where the error is the statistical one. 

Two systematic uncertainties may affect this analysis, 
the first related to the assumed primary CR composition. 
In Fig. [T3]the dependence of the Moon shadow displace- 
ment on the fraction of protons in the primary spectrum 
is shown as a function of the multiplicity. The proton ra- 
tio has been varied by ±10% from the assumed standard 
chemical composition. Indeed, in the investigated energy 
range, more than 90% of the CRs triggering ARGO-YBJ 
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FIG. 15: Expected westward displacement of the Moon 
shadow as a function of multiphcity, calculated assuming dif- 
ferent primary composition models. The dashed curves show 
the 7% shift, corresponding to (Jchem, from the solid line (see 
text). The multiplicity bins are shown by the horizontal bars. 
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FIG. 16: Expected westward displacement of the Moon 
shadow as a function of multiplicity calculated assuming dif- 
ferent hadronic interaction models. The solid and dashed 
curves are the best-fit results assuming the QGSJET (upward 
black tringles) and SIBYLL (downward red triangles) models. 
The multiplicity bins are shown by the horizontal bars. 



are protons and He nuclei [27|, whose spectra have been 
measured with uncertainties less than 10% [23|. The re- 
sults have been fitted with function (j31) obtaining the 
systematic uncertainty associated to the chemical com- 
position, CTc/iem = 7%. 

The second source of systematic uncertainty may be 
related to the use of different hadronic interaction mod- 
els. The results obtained with the QGSJet and SIBYLL 
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FIG. 17: Upper panel: Displacement of the Moon shadow 
from the apparent center in the North-South direction as a 
function of the observation time. Middle panel: Displacement 
of the Moon shadow from the apparent center in the East- 
West direction as a function of the observation time. Lower 
panel: The ratio of the observed deficit count to the expected 
one as a function of the observation time. The plots refer to 
events with a multiplicity Nstrip >100. 



codes are compared in Fig. [1^1 The different displace- 
ments have been fitted with function ([3]) obtaining the 
systematic uncertainty associated to these models, Ohadr 
— 12%. We expect that this uncertainty will be reduced 
by using new hadronic codes developed on the basis of 
the LHC data. 

Finally, the difference in the energy dependence of 
the Moon shadow displacement between data and Monte 
Carlo simulation has been estimated to be -1-5% ± 8stat% 

± 7c/iem% ± (12/ia(;r/2)%. 

The absolute rigidity scale uncertainty in the ARGO- 
YBJ experiment is estimated to be smaller than 

13% in the en- 
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ergy range from 1 to 30 (TeV/Z), where the Moon shadow 
is shifted from its position due to the effect of the GMF. 
For 7-induced showers we expect a lower scale uncer- 
tainty due to lack of uncertainties related to the hadronic 
interaction models and to the primary composition. 



D. Long-term stability of the detector 

The stability of the detector performance as far as 
the pointing accuracy and the angular resolution are 
concerned is crucial in 7-ray astronomy. Since Novem- 
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ber 2007 the full detector is in stable data taking with 
duty cycle >85%. Therefore, the stability of the ARGO- 
YBJ experiment has been checked by monitoring both 
the position of the Moon shadow, separately along R.A. 
and DEC. projections, and the amount of shadow deficit 
events in the period November 2007 - November 2010, 
for each sidereal month and for events with N strip >100. 

As discussed in Sec. V, B, the displacement of the 
center of the Moon shadow in the North-South direction 
enables us to estimate the systematic error in pointing 
accuracy and its long-term stability aside from Monte 
Carlo simulations, since the East- West component of the 
GMF is almost zero at YangBaJing. The displacement 
of the shadow position from the Moon center in the 
North- South direction is plotted in the upper panel of 
Fig. [T7] as a function of the observation time. Assuming 
a constant function, the best-fit result (continuous line) 
shows that the Moon shadow is shifted towards North by 
(0.19±0.02)°. The RMS around this position is 0.13°. 

In the middle plot the displacement along the East- 
West direction is shown. The best-fit result (continuous 
line) shows that the Moon shadow is shifted towards West 
by (-0.36±0.02)°, in agreement with the Monte Carlo ex- 
pectations ((-0.35±0.07)°). The RMS around this Moon 
position is 0.11°. 

The amount of CR deficit due to the Moon provides 
a good estimation of the size of the shadow, therefore of 
the angular resolution. 

The observed number of deficit events Ndef{< R) 
within an angular distance R from the Moon center is 
approximately given by 



Ndef{<R)= [l-e-°-^-(^ 



■N„ 



(4) 



where Nmoon is the number of events intercepted by the 
Moon and a is the Gaussian width of the shadow. N^oon 
is estimated by fitting the Eq. 2] to the experimental 
points fixing ct=0.8°. The expected deficit events are 
simply counted from the background content within the 
Moon disc. In the lower plot of Fig. [T7] the ratio of the 
observed deficit count to the expected one is shown as a 
function of the observation time. The line shows the best- 
fit result assuming a constant function: 1.005±0.016. 
The RMS of the corresponding distribution is 0.11. 

In conclusion, the position of the Moon shadow mea- 
sured with the ARGO-YBJ experiment turned out to be 
stable at level of 0.1° and the angular resolution is stable 
at level of 10%, on a monthly basis. 



E. The day/night effect 

Some differences in the Moon shadow signal could arise 
due to the interaction of the GMF with the solar wind. 
We can point them out by dividing the data sample ac- 
cording to the relative positions of the Moon-Sun system, 
i.e. by requiring that the angular distance between the 
Moon and the Sun is smaller or larger than 90°. This 



selection can be carried out by dividing all recorded data 
in two "day" /"night" samples. 

A similar analysis has been performed by the Ti- 
bet AS7 array using 1990 - 1992 data ^3^ and by the 
MACRO underground detector with high-energy muons 
in the period from 1989 to 2000 [3], though at higher ener- 
gies (Tibet AS7 at -10 TeV and MACRO at -20 TeV) 
and with less statistics. While the Tibet collaboration 
did not find any day-night effect, the MACRO experi- 
ment observed a sharper shadow in the "night" sample 
and a broader shadow in the " day" one, with a lower sig- 
nificance, concluding that the night events encounter a 
reduced magnetic field with respect to the day events. 

The results of the analysis carried out with the ARGO- 
YBJ experiment for the two subsamples are shown in 
Fig. [TH] for events with Nstrip >60. The data refer to 
the November 2007 - November 2009 period and in each 
sample the Moon is observed for about 1150 hours. We 
did not find any appreciable difference between the "day" 
and the "night" shadows (statistical significance of the 
maximum deficit 31 vs 30 s.d.) shadows. Accordingly, 
the shape of the Moon shadow seems to be independent 
of the position of the Moon with respect to the Sun. This 
implies that effects due to the solar wind do not give a 
considerable contribution to the CR bending, at least in 
the period of minimum of the solar activity. 



VI. CONCLUSIONS 

The Moon shadowing effect on cosmic rays has been 
observed by the ARGO-YBJ experiment in the multi- 
TeV energy region with a statistical significance greater 
than 55 standard deviations. 

We observed a westward displacement, due to the 
GMF, up to about 1.3°, proving the detection of the 
Moon shadow cast also by sub- TeV primaries. By means 
of an accurate Monte Carlo simulation of the CR propa- 
gation in the Earth-Moon system we have studied in de- 
tail the role of the GMF and of the detector PSF on the 
observed shadow. The measured deficit counts around 
the Moon position are found in fair agreement with the 
expectations based on the primary cosmic ray composi- 
tion derived from the direct observational data. 

The dependence of the measured angular resolution on 
the particle multiplicity is in good agreement with Monte 
Carlo calculations. A systematic shift of (0.19±0.02)° 
towards North has been observed. 

We have estimated the primary energy of the detected 
showers by measuring the westward displacement as a 
function of the multiplicity, thus calibrating the relation 
between shower size and CR energy. The systematic un- 
certainty in the absolute rigidity scale is evaluated to be 
less than 13% in the range from 1 to 30 (TeV/Z), mainly 
due to the statistical one. 

The position of the Moon shadow measured with the 
ARGO-YBJ experiment turned out to be stable at a level 
of 0.1° and the angular resolution stable at a level of 10%, 
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on a monthly basis. These results make us confident 
about the detector stability in the long-term observation 
of gamma-ray sources. 

Finally, we have studied with high statistical accuracy 
the shadowing effect in the " day" /" night" time looking 
for possible effects induced by the solar wind. Within the 
statistical accuracy of this study we find that the solar 
wind does not give appreciable contribution to the CR 
bending, at least in the period of minimum of the solar 
activity. 
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Appendix 

In this Appendix the analytical calculation of Eq. ([1]) 
is presented. Since only the magnetic field is supposed 
to act upon the particles trajectories, they read as: 

Zec^ /■* r , dx _, , , , 

x(t) = xq -I- vo t -I -— / dr / da-— xB(x, a) (5) 

E Jq Jo da 

at time i in a certain reference frame, where: 

• x(t) is the particle position at time t; 

• xq and vq are the initial position and velocity of 
the particle; 

• Ze and E are its charge and its (constant) energy; 

• B(x, t) is the magnetic field, which in principle can 
vary with respect to both position and time; 

• a is the time the inner integral is computed over. 

If it is possible to write an explicit functional form for 
B(x, i), an attempt to solve Eq. ([5]) can be made. On 
the contrary, especially when no analytical expression of 
the time behavior is known, the equation can be solved 
with numerical techniques. 



Equation ([5]) explicitly shows the perturbation induced 
by the magnetic field on the straight trajectory (x(i) = 
Xq + vgi) . It suggests an iterative method to determine 
the solution, which can be expressed as the series: 

x(i) =xo(B«){t) +xo(Bi)(i) + ... 
where xo(BO)it) = xq -I- vpt is the unperturbed 
(straight) trajectory and for higher orders we have: 



Ax, 



Zee' 



dx. 



O(B') 



da 



O(S'+i)(0 - 
X B(xo(i3i),Q) 

i = 0,l,... 



where AxQ(gi+i^{t) — XQ(gi+i^{t) — (xq + vq t) is the 
displacement from the unperturbed trajectory at time t. 
At the first-order approximation we find: 

Zee' /■* r 

Ax(i) ~ — p;— vo X / dr / daB(xo +Voa,a) 
E Jo Jo 



Ax(i) ~ — vo xTB(i;xo,vo) 



(6) 



where lB(i;xo,vo) is the integral of the magnetic field 
along the straight trajectory, whose value depends only 
on the time of the motion (t) and on its initial conditions 
(xo and Vq). 

Since the phenomenon studied concerns ultrarelativis- 
tic particles, once we fix the initial position and the final 
time, Eq. © reads: 



Ax: 



— Vo X 2'b(vo) 



At the first approximation the displacement depends only 
on the charge-to-energy ratio of the primary and on the 
initial direction of its ultrarelativistic motion (versor vq). 
Now, let us consider only the lowest order of the GMF 
multipoles-expansion, i.e. the dipole term: 



B(x) 



3(b • x)x - x'h 



where b has intensity b w 8.1 • 10^"^ T m'^ and the south 
magnetic pole is supposed to have coordinates 78.3° 
South, 111.0° East. By setting vo||xo {vertieal direetion 
approximation) and integrating from YangBaJing to a 
distance ^ 60 Earth radii, Eq. ([1]) is immediately ob- 
tained. 
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FIG. 18: Significance map of the Moon region observed during 
"day" (panel (a)) and "night" (panel (b)) for events with 
Nstrip >60. The color scale gives the statistical significance 
in terms of standard deviations. 



